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Abstract 

This paper is the first application of the compensation approach (a well-established the- 
ory in probability theory) to counting problems. We discuss how this method can be ap- 
plied to a general class of walks in the quarter plane 1? + with a step set that is a subset of 
{(—1, 1), (—1, 0), (—1, —1), (0,-1), (1, —1)} in the interior of We derive an explicit expres- 
sion for the generating function which turns out to be nonholonomic, and which can be used to 
obtain exact and asymptotic expressions for the counting numbers. 
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1 Introduction 

In the field of enumerative combinatorics, counting walks on a lattice is among the most classical 
topics. While counting problems have been largely resolved for unrestricted walks on Z 2 , walks 
that are confined to the quarter plane still pose considerable challenges. In recent years, much 
progress has been made, in particular for walks in the quarter plane with small steps, which means 
that the step set S is a subset of : \j\ ^ 1} \ {(0,0)}. Bousquet-Melou and Mishna [7] 

constructed a thorough classification of these models. By definition, there are 2 8 such models, 
but after eliminating trivial cases and exploiting equivalences, it is shown in [7] that there are 79 
inherently different problems that need to be studied. Let qi & denote the number of paths in 
having length k, starting from (0,0) and ending in and define the generating function (GF) 

as 

oo 

Q(x,y;z)= ^ q iJtk x l y J z k . (1) 

i ,j,k=0 

There are then two key challenges: 

(i) Finding an explicit expression for Q(x,y;z). 
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(ii) Determining the nature of Q(x,y;z): is it holonomic (the vector space over C(x,y,z) — the 
field of rational functions in the three variables x, y, z — spanned by the set of all derivatives 
of Q(x,y;z) is finite dimensional, see [11, Appendix B.4])? And in that event, is it algebraic, 
or even rational? 

The common approach to address these challenges is to start from a functional equation for the 
GF, which for the walks with small steps takes the form (see [7, Section 4]) 

K(x, y; z)Q(x, y; z) = A{x)Q(x, 0; z) + B(y)Q(0, y; z) - 5Q(0, 0; z) - xy/z, (2) 

where K, A and B are polynomials of degree two in x and/or y, and 5 is a constant. For z = 1/|<S|, 
(2) belongs to the generic class of functional equations (arising in the probabilistic context of 
random walks) studied and solved in the book [8]. For general values of z, the analysis of (2) for 
the 79 above-mentioned models has been carried out in [13, 18], which settled Challenge (i). 

In order to describe the results regarding Challenge (ii), it is worth to define the group of 
the walk, a notion introduced by Malyshev [15]. This is the group of birational transformations 
W = (£,77), with 

(M = LgBcM±\, = (lpW, A (3 ) 

V y l^(i,+i)es x J \ x l^(+i,j)esy J J 

which leaves invariant the function Yl(ij)eS X% V 3 ■ Clearly, £ o £ = 77 o rj = id, and TV is a dihedral 
group of even order larger than or equal to four. 

Challenge (ii) is now resolved for all 79 problems. It was first solved for the 23 models that 
have a finite group. The nature of the GF (as a function of the three variables x, y, z) was de- 
termined in [7] for 22 of these 23 models: 19 models turn out to have a GF that is holonomic 
but nonalgebraic, while 3 models have a GF that is algebraic. For the 23rd model, defined by 
S = {(—1, 0), (—1, —1), (1, 0), (1, 1)} and known as Gessel's walk, it was proven in [5] that Q(x, y; z) 
is algebraic. Alternative proofs for the nature of the GF for these 23 problems were given in [9]. For 
the remaining 56 models, which all have an infinite group, it was first shown that 5 of them (called 
singular) have a nonholonomic GF, see [16, 17]. Bousquet-Melou and Mishna [7] have conjectured 
that the 51 other models also have a nonholonomic GF. Partial evidence was provided in [6, 18], 
and the proof was recently given in [14]. 

In this paper we consider walks on with small steps that do not fall into the class considered 
in [7]. The classification in [7] builds on the assumption that the steps on the boundaries are the 
same steps (if possible) as those in the interior. However, when the behavior on the boundaries is 
allowed to be different, we have many more models to consider. There are several motivations to 
permit arbitrary jumps on the boundary. First, this allows us to extend the existing results, and 
thus to have a better understanding of the combinatorics of the walks confined to a quarter plane. 
Second, some studies in probability theory (see e.g. [12]) suggest that the behavior of the walk on 
the boundary has a deep and interesting influence on the quantities of interest. In this paper we 
first consider the walk in Figure 1, with steps taken from S = {(—1,1), ( — 1,-1), (1,-1)} in the 
interior, Sjj = {(— 1, 1), (— 1, 0), (1, 0)} on the horizontal boundary, Sy = {(0, 1), (0, — 1), (1, — 1)} 
on the vertical boundary, and Sq = {(0, 1), (1,0)} in state (0,0). 

Notice that the step set S = {(— 1,1),(— 1,-1), (1,-1)} in the framework of [7] would render 
a trivial walk, since the walk could never depart from state (0,0). However, by choosing the steps 
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Figure 1: The jumps of the walk 

on the boundaries as in Sh, Sv an d Sq, it becomes possible to start walking from state (0,0), and 
we have a rather intricate counting problem on our hands. 

It turns out that our walk has an infinite group. To see this, observe that the interior step set 
in Figure 1 and the one represented in Figure 2 have isomorphic groups ([7, Lemma 2] says that 
two step sets differing by one of the eight symmetries of the square have isomorphic groups), and 
notice that the group associated with the step set in Figure 2 is infinite ([7, Section 3.1]). 
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Figure 2: Walk considered in [17] 

Since our walk has an infinite group, the approach in [7] cannot be applied. The methods 
developed in [17] also fail to work. Indeed, the main tool used there is an expression of Q(x, 0; z) and 
Q(0, y; z) as series involving the iterates of the roots of the kernel. While these series are convergent 
when the transitions (—1,0), (—1,-1), (0,-1) are absent, they become strongly divergent in our 
case, see [8, Chapter 6]. The approach via boundary value problems of [8] seems to apply, but this 
is more cumbersome than in [13, 18]. Indeed, contrary to the 79 models studied there, for which 
in (2), A and B depend on one variable, and 5 is constant, these quantities are now polynomials 
in two variables, since the functional equation (2) for our walk becomes (for a proof, see Equations 
(12)-(17) in Section 2) 

K(x, y; z)Q(x, y; z) = [l+x 2 -x 2 y-y]Q(x, 0; z) + [l+y 2 -xy 2 -x]Q(0, y; z)+[x+y-l]Q(0, 0; z)-xy/z 

(4) 

with the kernel 

K(x, y; z) = 1 + x 2 + y 2 — xy/z. (5) 

We shall derive an explicit expression for Q(x,y;z), which turns out to be a meromorphic 
function of x, y with infinitely many poles in x and in y. This implies that Q(x, y; z) is nonholonomic 
as a function of x, and as a function of y. We note that for one variable, a function is holonomic if 
and only if it satisfies a linear differential equation with polynomial coefficients (see [11, Appendix 
B.4]), so that a holonomic function must have finitely many poles, since the latter are found among 
the zeros of the polynomial coefficients of the underlying differential equation. Accordingly, the 
trivariate function Q(x, y; z) is nonholonomic as well, since the holonomy is stable by specialization 
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of a variable (see [11]). As in [14, 18], this paper therefore illustrates the intimate relation between 
the infinite group case and nonholonomy. 

The technique we are using is the so-called compensation approach. This technique has been 
developed in a series of papers [1, 3, 4] in the probabilistic context of random walks; see [2] for an 
overview. It does not aim directly at obtaining a solution for the GF, but rather tries to find a 
solution for its coefficients 

oo 

%j( Z ) = ^2<lij,k zk - (6) 
fc=0 

These coefficients satisfy certain recursion relations, which differ depending on whether the state 
lies on the boundary or not. The idea is then to express qij(z) as a linear combination of 
products a*/3 J , for pairs (a, (5) such that 

K{\fa,\/foz)=Q. (7) 

By choosing only pairs (a,/3) for which (7) is satisfied, the recursion relations for qij(z) in the 
interior of the quarter plane are satisfied by any linear combination of products a l /3 J , by virtue of 
the linearity of the recursion relations. The products have to be chosen such that the recursion 
relations on the boundaries are satisfied as well. As it turns out, this can be done by alternatingly 
compensating for the errors on the two boundaries, which eventually leads to an infinite series of 
products. 

This paper is organized as follows. In Section 2 we obtain an explicit expression for the gen- 
erating function Q(x, y; z) by applying the compensation approach. In Section 3 we derive an 
asymptotic expression for the coefficients qij t k f° r large values of k, using the technique of singular- 
ity analysis. Because this paper is the first application of the compensation approach to counting 
problems, we also shortly discuss in Section 4 for which classes of walks this compensation approach 
might work. 



2 The compensation approach 

We start from the following recursion relations: 

%j,k+l = 9i-ij'+l,fc + Qi+l,j-i,k + Qi+l,j+l,ki i,j^l,k^ 0, (8) 

9i,0,fc+i = + <Zi+l,l,fc + Qi-l,0,k + <?i+l,o,fc; i ^ l,k ^ 0, (9) 

Qo,j,k+i = Qi,j-i,k + ?i,j+l,fc + Qo,j-l,k + 9o,i+i,fe, 3 5 s 1, k > 0, (10) 

go,o,fc+i = ?o,i,fc + gi.i.fc + gi.o.fc' k ^ o. (11) 

Since qij,o = if i + j > and (/o,o,o = 1> these relations uniquely determine all the counting 

numbers qijk- Multiplying the relations (8)-(ll) by z k and summing w.r.t. k ^ leads to (with 
the generating functions q^j = qij(z) defined in (6)) 

Qi,j/z = Qi-iJ+i + + i,3 > 1, (12) 

qi,o/z = gi-i,i + qi + i,i + gi-i,o + gi+i,o, i > 2, (13) 

qo,j/z = gij-i + gij + i + qoj-i + goj+ii 3 > 2 , (14) 
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and 



?i,oA = 9o,i + 92,1 + 9o,o + 92,0, (15) 
9o,i/2 = qifi + ^1,2 + 9o,o + 9o,2, (16) 
qo,o/z = l/z + g ,i + 9i,i + 9i,0- (17) 

(For a reason that will be clear below, we want to discus separately the equations involving qo,o, 
and we therefore do not merge (15) (resp. (16)) and (13) (resp. (14)). 

Lemma 1. Equations (12) -(17) have a unique solution in the form of formal power series. 

Proof The generating function qij can be written as a Taylor series about z = 0. Substituting 
these Taylor series into Equations (12)-(17), and equating coefficients of z k , yields Equations (8)- 

(11) . The latter obviously have a unique solution for the coefficients 9ij,fc, because these counting 
numbers can be determined recursively using go, 0,0 = 1- d 

In order to find the unique solution for qtj, we shall employ the compensation approach, which 
consists of three steps: 

• Characterize all products a l (3 3 for which the inner equations (12) are satisfied, and construct 
linear combinations of these products, which in addition to being formal solutions to (12), 
also satisfy (13) and (14). 

• Prove that these solutions are formal power series. 

• Determine the complete unique solution qi j by taking into account the boundary conditions 
(15)-(17). 

2.1 Linear combinations of products 

Substituting the product a l j3 3 into the inner equations (12), and dividing by common powers, yields 

a/3/z = p 2 + a 2 + a 2 p 2 . (18) 

Incidentally, note that (18) is nothing else but (7) (see (5)). Hence, a product a i P J is a solution of 

(12) if and only if (18) is satisfied, and any linear combination of such products will satisfy (12). 
Figure 3 depicts the curve (18) in for z = 1/4. 

Now we construct a linear combination of the products introduced above, which will give a 
formal solution to the equations (12)-(14). The first term of this combination, say olqPq, we require 
to satisfy both (12) and (13). In other words, the pair (ao. A)) nas to satisfy (18) as well as 

a p/ z = p 2 + a 2 f3 2 + + a 2 13. (19) 

The motivation to start with a term satisfying both (18) and (19) will be explained at the end of 
this subsection, see Remark 5. 

Lemma 2. There exists a unique pair ((Xq,Pq), with 



l-VT^s^ _ a 2 
«o — — , Po — — — 2' v 2U J 



4z l + ag 
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Figure 3: The curve (18) in R 2 + for z = 1/4. 



of formal power series in z such that the product ckq/Sq satisfies both (12) and (13). 
Proof. Equation (19) yields 

a/z = (3 + a 2 f3 + 1 + a 2 = (1 + /3)(1 + a 2 ). (21) 

Subtracting (19) from (18) gives /3(1 + a 2 ) = a 2 , so that (1 + /3)(1 + a 2 ) = 1 + 2a 2 , and together 
with (21) this gives a/z = 1 + 2a 2 . The last equation has the two solutions (1 ± yl — 8z 2 )/(4z), 
and only the solution with the — sign is a formal power series in z. □ 

The function Oq/Sq thus satisfies (12) and (13), but it fails to satisfy (14). Indeed, if a l (3 3 is a 
solution to (14), then (a,/3) should satisfy 

a/3/z = a 2 + a 2 /3 2 + a + a/3 2 , (22) 

and (ao,/3o) is certainly not a solution to (22). 

Now we start adding compensation terms. We consider coa /3Q + d\a % ^ d , where (a,/3) satisfies 
(18) and is such that aja^o + ^ia*/3 J satisfies (14). The identity (14), which has to be true for any 
j ^ 2, forces us to take (3 = (3q. Then, thanks to (18), we see that a is the conjugate root of ao i n 
(18). Denote this new root by a±. Hence, ao and a\ are the two roots of Equation (18), where /3 
is replaced by /3o- In other words, ao and a\ satisfy 

(1 + f3 2 )a 2 - {Po/z)a + /3 2 = 0. (23) 

Due to the root-coefficient relationships, we obtain 

^L = ao + ai . (24) 

We now determine d\ in terms of cq. Note that coa /3Q + d\a\P 3 Q is a solution to (14) if and only if 
(A)/z)(c + di) = (a c + aidi)(l + Pi) + (cq + di)(l + /3 2 ). (25) 
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The latter identity can be rewritten as 



f° Z R2 (co + di) = (a c + a\di) + (c + d±), (26) 

and thus, using (24), 

di = ^cq. (27) 
1 - a 

The series coa pQ + ^i a iPo & fter one compensation step satisfies (12) and (14) for the interior and 
the vertical boundary. However, the compensation term d\a\f3Q has generated a new error at the 
horizontal boundary. To compensate for this, we must add another compensation term, and so on. 
In this way, the compensation approach can be continued, which eventually leads to 

H H H 



Xi,j = c a Po +d 1 a\P J + c\a\P>{ + d 2 a l 2 fi{ + c 2 a 2 /3 J 2 H (28) 



v v 
where, by construction and (18), for all k ^ we have 



1 _ /T _ A z 2(l _j_ +2) 

Pk = f(<* k ), a k+1 = f(f3 k ), f(t) = V 2z{1 + t 2) } t, (29) 

and the function f{t) follows from (18). The construction is such that each term in (28) satisfies 
the equations in the interior of the state space, the sum of two terms with the same a satisfies 
the horizontal boundary conditions (H) and the sum of two terms with the same f3 satisfies the 
vertical boundary conditions (V). Hence Xij satisfies all equations (12)-(14). See also Figure 4. 
Furthermore, 

1 - a k+ l 1 - Pk+l j , \ n (nn\ 

d k+ i = c fc , c k+1 = —d k+1 , k^O. (30) 

1 - OL k 1 - Pfe 

An easy calculation starting from (30) yields 

(l-Qfc+lXl-pfc) (1 ~ gfc+lKl - Qfc+i) , . n 

4+1 = 7z rr- FT-cq, c fc+ i = — — — — c , k ^ 0, (31) 

(1 - a )(l - Po) (1 - «o)(l - Po) 

so that choosing (arbitrarily) cq = (1 — ao)(l — Po) finally gives 

oo 

x ij = I> - WK 1 - «fc)4 - (1 - «fc+i)4+i]- (32) 

fc=0 

By symmetry, we also obtain that formal solution to (12)-(14). This leaves to consider 

the solution for the states (0,1), (0,0) and (1,0), i.e., the boundary conditions (15)-(17). 
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Figure 4: The curve (18) near for z = 1/4, and the construction of the products. 
2.2 Formal power series 

Proposition 3. The sequences {ctfclfc^o and {/3fc}fc^o appearing in (28) -(29) satisfy the following 
property: for all k ^ ; a& = z 2k+1 cik(z) and ftk = z 2k+2 f3k{z) , where and /3f. are formal power 
series such that 2^(0) = /3fc(0) = 1. 

Proof. For / defined in (29), all p ^ 0, and all real numbers s±, s%, . . ., 

f(z p [l + 8l z + s 2 z 2 + •••]) = z p+1 [l + siz + (1 + s 2 )z 2 + •••]. (33) 

Because (20) yields ao = z + 2z 3 + • • • , the proof is completed via (29). □ 

The following result is an immediate consequence of Proposition 3. 

Corollary 4. For i,j ^ 0, Xij and Xj t i defined in (32) are formal power series in z. 

Remark 5. The approach outlined in Subsection 2.1 is initialized with a term satisfying both (18) 
and (19). Alternatively, we could start with an arbitrarily chosen term with power series «o an d 
/3o satisfying (18) only. This term would violate (19) as well as (22), and therefore generate two 
sequences of terms, one starting with compensation of (19) and the other with (22). It is readily 
seen that in one of the two sequences, some of the parameters a or /3 will not be power series. 
Indeed, the two sequences of terms would be generated as in (29), with instead of / the function 



f (f) - i±vi-4* a (i+o 

f±{t) ~ WT^) *■ 



(34) 



In Subsection 2.1, we already described the construction of the sequence with the choice /_ = /. 
For / = /+, notice that for all p ^ 0, and all real numbers si, . . ., we have 



f + ( z P[l + s 1 z + ---]) = zP~ 1 [l + s 1 z + ---[ 



(35) 



Hence, the resulting solution constructed as in (29) would fail to be a formal power series, cf. Equa- 
tion (32) and Subsection 2.3. 
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2.3 Determining the unique solution 

We now determine the generating functions qij. Define 

%ij — (36) 

so that in particular 



xo,o = 2 X/( 1 _ Pk)(<*k+i ~ a k)- (37) 

k=0 



Proposition 6. The expressions 

and 
with 



9o,o = z , 7z (38) 

l-2z + ZX n 



qij=cxij, i + j>0, (39) 
1 



(40) 



1 — 2z + zxo t o ' 

are the unique solutions to Equations (12)-(17) in the form of formal power series. 
Proof. We shall look for a solution satisfying, for i + j > 0, 

By symmetry we must have c = c. Such combinations have been shown to satisfy (12)-(14), 
since they do not involve go,o- We should still determine the solution for which also the boundary 
conditions (15) and (17) (or by symmetry (16)) are satisfied. This can be achieved by choosing c 
and (/o.o appropriately. 

By using (17) and (15) (or (17) and (16)), we obtain that 

c = - pz - zA, oT^ - ( 42 ) 

£1,0 - Z[Xo,l + X2,l + X 2 fl\ ~ Z z [Xlfl + X ,l + Xl,lJ 

and 

9o,o = 1 + zc[x 0j i + x lt x + x lj0 ]. (43) 
In addition, tedious calculations starting from (32) give 

xi >0 /z - [x ,i + x 2 ,i + x 2 ,o + x ,o] = 1, (V 2 - l)a?o,o + 2 - [xi, + x ,i + £1,1] = 0. (44) 

Thanks to the equations in (44), the denominator of (42) can be written as z[l — 2z + .zSo.o]- The 
functions defined in (38) and (39) are the solutions to Equations (12)-(17) by construction. Finally, 
the uniqueness follows from Lemma 1. □ 
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Proposition 7. 

oo 

Q{x,y;z) = ^ qijx'y 3 (45) 

1 - a k 1 - a k+1 



k=0 



1-/3, 



k 



PkV 



1 - a k x 1 - afc+ixj 1 - (3 k x 



1 - a k y 1 - afc+iy 



(46) 



with qij as in Proposition 6, c as in (40), and a k and (3 k as in (20) and (29), respectively. Fur- 
thermore, 

Q(0,0;z) =c(l + x ,o), Q(l,0;2t)=Q(0,l^) = c(l-oo(«)), Q(l,l;*)=c (47) 

Proof. This is immediate from (32), (37)-(40). □ 

Proposition 8. For all < \z\ < 1/3, i/ie functions YlTj=o Qi,j xl y J > Yl'tLo Qi,o x% an( ^ YlJLo Qo,jy J 
have infinitely many poles. 

Proof. This is a direct consequence of Propositions 3 and 7. □ 

It follows immediately from Proposition 8 that both J^£o 1ifi xt ano - Yl'jLo lojU' 1 are nonholo- 
nomic. The holonomy being maintained after specialization to a variable (see [11]), we reach the 
conclusion that the trivariate generating function Q(x,y;z) is nonholonomic as well. 

2.4 Retrieving coefficients 

Now that we have an explicit expression for q$j, we briefly present an efficient procedure for 
calculating its coefficients Qijk- 

To compute qij we need to calculate in principle an infinite series, or to solve the recurrence 
relations. However, we shall prove that if we are only interested in a finite number of coefficients 
Qij ki then it is enough to take into account a finite number of a k and f3 k . For k ^ 0, denote by 
k V 1 the maximum of k and 1. Define 

N^ = l + L±max{p-(iVl + 2(jVl)),p-(2(*Vl)+jVl)}J . (48) 

Proposition 9. For any i,j ^ 0, the first p coefficients of qij only require the series expansions 
of order p of a , A), • • ■ , fi^ . 

Proof. In order to obtain the series expansion of qo,o of order p, it is enough to know the series 
expansion of xo,o of order p, see (38). From (32) we obtain 

oo oo 

sb,o = 2^(1 - /9fc)(ttfc+i - a k ) = -2q + 2^/3 fc (a fc - a k+l ). (49) 

k=0 k=0 

With Proposition 3, f3 k a k = 0(z 4fc+3 ) and f3 k a k+ i = 0(z Ak+5 ), so that in order to obtain the series 
expansion of xo,o of order p, it is enough to consider in (49) the values of k such that Ak + 3 ^ p. In 
other words, we have to deal with a®, (3q, . . . , a k , f3 k , a k+ \ for 4/c+3 ^ p. Since Np'° = 1+ [(p— 3)/4j , 
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see (48), Proposition 9 is shown for i = j = 0. The proof for other values of i and j is similar and 
therefore omitted. □ 



As an application of Proposition 9, let us find the numbers (?o,o,fc for k € {0, . . . , 10}. Taking i = 

j = 0, we have A{q° = 2; we then calculate the series expansions of order 10 of ao, Po> Pi,ot2, 02- 
After using (38) and (49), we obtain the numbers: 

1, 0, 2, 2, 10, 16, 64, 126, 454, 1004, 3404. (50) 



3 Asymptotic analysis 

In this section we derive an asymptotic expression for the coefficients qij,k f° r large values of k, using 
the technique of singularity analysis (see Flajolet and Sedgewick [11] for an elaborate exposition). 
This requires the investigation of the function qij = qij(z) near its dominant singularity (closest 
to the origin) in the 2-plane. First, notice that 

oo 
i,j=0 

because there are at most 3 fc paths of length k. This clearly implies that qij,k ^ 3 fc , so that qij is 
analytic at least for \z\ < 1/3, and the singularities z of qij must satisfy \z\ ^ 1/3. 

In fact, the singularities of qij are given by the singularities of Xij and the zeros of the denom- 
inator of c in (40), see (38) and (39) . In what follows, we first show in Subsection 3.1 that x-ij 
converges in the domain \z\ ^ l/\/8- We then find in Subsection 3.2 the first singularity of qij. 
Finally, we state and prove our asymptotic results in Subsection 3.3. 



3.1 Convergence of the solutions 

Since Xij is constructed from the functions aoj A)j a i> ■ ■ ■■> an d since all these functions follow 
from the iterative scheme (29), i.e., from ao, it is clear that the singularities of Xij are in fact the 
singularities of ao, namely (see (20)) 

Z = ± -L. (52) 

Further, we have that: 

Proposition 10. The sequences {a^lfc^o an d {Pk}k^o appearing in (28)-(29) satisfy the following 
properties: 

(i) For all \z\ < 1/V8, we have l/y/2 ^ |a | > |/?o| > |"i| > > ••• ■ 

2fc+l 2k-\-2 

(ii) For all \z\ l/v8 and all k 0, we have |a^| ^ l/v2 and \/3k\ ^ l/v2 

Before starting with the proof of Proposition 10, let us recall Rouche's theorem (see e.g. [11, 
Chapter 4]). 

Theorem 11 (Rouche's theorem). Let the functions F and G be analytic in a simply connected 
domain of C containing in its interior the closed simple curve 7 . Assume that F and G satisfy 
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\F(r) — G(r)\ < \G(r)\ for r on the curve 7 . Then F and G have the same number of zeros inside 
the interior domain delimited by 7 . 

Proof of Proposition 10. To show (i) and (ii), it is enough, thanks to (29) and since ao = to 
prove that for all \t\ ^ 1 and \z\ ^ l/\/8j 

\f(t)\ < J|. (53) 

In order to show (53), we first note that / satisfies the algebraic equation (l+£ 2 )/(t) 2 — tf{t)/z+t 2 = 
0, see (18) and (29). As a consequence, the function f(t)/t is such that (with r = f{t)/t) 

{l+t 2 )r 2 -r/z + 1 = 0. (54) 

Let us now fix \t\ ^ 1, \z\ < l/\/8, and introduce the functions 

F(r) = (1 +t 2 )r 2 - r/z + 1, G(r) = -r/z + 1. (55) 

We have 

\F(r)-G(r)\ = \(l + t 2 )r 2 \^2\r\ 2 , \G(r)\ = \ - r / z + 1\ > \r\/\z\ - 1. (56) 
Further, a straightforward computation gives that for any 

1 "^T W <'-»< 1 + V ^T W («) 

4|z| 4|z| 

we have, for all \r\ = r z , 

2\r\ 2 < \r\/\z\ - 1. (58) 

It then follows from Rouche's theorem, applied to the circle with radius r z , that (54) has one unique 
solution r = f{t)/t with \r\ < r z . By a continuity argument, we deduce that for all \t\ ^ 1 and all 

\z\ < l/y/8, 

m < iz^EM. (59) 

Equation (53) and Proposition 10 follow immediately. □ 

Corollary 12. For i,j ^ 0, Xij and xj^ defined in (32) are convergent for \z\ ^ l/\/8. Further, 
for each \z\ ^ l/\/8, 

00 

J] 1x^(^)1 < 00. (60) 
i,i=o 

Proof. The first statement in Corollary 12 immediately follows from Proposition 10 and from the 
expression for Xij. Next, use (32) to get that for i + j > 0, 



00 

fc=0 



j\ ^J^]! - pk\\Pk\ j [\l-akWakt + \l-a k+l \\a k+1 t] . (61) 
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Using then that |1 — /3k\, |1 — aik\ an d |1 — afc+i| are bounded, see Proposition 10 (i), and applying 
Proposition 10(h) gives that there is some constant constant C > such that for all i + j > 0, 



\ Xi>j \^Cj2l/V2 {2k+1)i+{2k+2)j . (62) 

k=0 

Likewise, we can prove that for some constant C > 0, 

oo 

h,„KC^l/^ +1) (63) 

fc=0 

The second claim in Corollary 12 immediately follows. □ 



3.2 Dominant singularity of q i; j 

A consequence of the results in Subsection 3.1 is that qij is analytic for \z\ < 1/V&, except possibly 
at points z such that the denominator of c in (40) equals 0. Denote this denominator by 

h(z) = l-2z + zxo, , (64) 

and let p denote the dominant singularity of c. 

Lemma 13. The series h(z) has a unique root in the interval (1/3, l/y/8), which is the dominant 
singularity p of c. 

The proof of Lemma 13 is presented in Appendix A. Furthermore, in Corollary 22 we shall 
prove that p G [0.34499975, 0.34499976] 

Remark 14. It seems hard to find a closed-form expression for p or to decide whether p is algebraic 
or not. From this point of view, the situation is quite different from that for the 74 (nonsingular) 
models of walks discussed in Section 1. Indeed, for all these 74 models, the analogue of p (which 
eventually will be the first positive singularity of Q(0, 0; z), Q(l, 0; z), Q(0, 1; z) and Q(l, 1; z), see 
below) is always algebraic (of some degree between 1 and 7), see [10]. 

3.3 Asymptotic results 

Here are our main results on the asymptotic behavior of large counting numbers: 
Proposition 15. The asymptotics of the numbers of walks qij,k as k —> oo is 

qi,j,k ~ C itj p~ k , (65) 

with p as defined in Lemma 13, 

c °>° = ' (66) 

and 
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Proof. Consider first the case i = j = 0. Thanks to (38), (64) and Lemma 13, we obtain that go,o 
has a pole at p, and is holomorphic within the domain 



{z € C : |*| < (1 + e)p} \ [p, (1 + e)p) (68) 

for any e > small enough. Moreover, the pole of (70,0 at p is of order one, as we show now. For 
this it is sufficient to prove that 1 + 5?o,o(p) an d that h'(p) 7^ 0, see again (38) and (64). In 
this respect, note that with (64) and Lemma 13 we have h(p) = 1 — 2p + pxo t o(p) = 0, so that 
1 + XQfi(p) = (3/9 — l)/p, which is positive by Lemma 13. On the other hand, it is a consequence 
of Lemma 23 that h'(p) 7^ 0. In particular, the behavior of (70,0 near p is given by 

1 + X 0fi (fi) 3/9-1 r . 

The holomorphy of go,o within the domain (68) and the behavior (69) of go,o near p immediately 
give the asymptotics (65) for i = j = (see e.g. [11]). The proof for other values of i and j is 
similar and therefore omitted. □ 

A result similar to Proposition 15 also holds for the total number of walks of length k (not 
necessarily ending in state (0,0)). Indeed, from formula (47), we obtain that 

Q(l,0;z)= Q(l,l;*) = -£r, (70) 

with h(z) as in (64). Accordingly, the exponential growth rate of these generating functions is the 
same as that of Q(0, 0; z), namely p (defined in Lemma 13). An identical reasoning as in the proof 
of Proposition 15 then gives the following result. 

Proposition 16. The asymptotics of the total number of walks YlTj=o1i,j,k as k — > 00 is 

00 1 

£ ^ ~ ^pm p " (71) 

i,j=0 y y ' ' 

The asymptotics of the total number of walks ending on the horizontal axis ( or, equivalently, on the 
vertical axis) Yli^o 1i,o,k as k ~ * 00 ^ 

V> 1 ~ «o(p) -k 

S^-^w' ' (72> 

Remark 17. A priori, the fact that 

g(0,0;z), Q(l,0;z), Q(0,l;z), Q(l,l;z) (73) 

have the same first positive singularity (namely, p) is far from obvious. This appears to be related 
to the fact that the two coordinates of the drift vector of our S are negative (they both equal —2). 
Here, for any step set S, the drift vector is given by (X]-i<i<i ^fj> S-i<j<i ^fj)> where Sfj = 1 if 
(i,j) € S, and otherwise. Indeed, the following phenomenon has been observed in [10] for the 
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79 models of walks mentioned in Section 1: if the two coordinates of the drift vector are negative, 
then the four functions (73) do have the same first positive singularity. 

Let us illustrate Proposition 15 with the following table, obtained by approximating p by 
0.34499975 and C 0i o by 0.0531, see Corollary 22 and Lemma 23 in Appendix A. 



Value of k 


Exact value of <7o,o,fc 


Approximation of qo,o,k 


Ratio 




(see Section 2.4) 


(see Proposition 15) 




10 


3.404 ■ 10 3 


2.222 • 10 3 


0.653 


20 


1.106 • 10 8 


9.305 • 10 7 


0.840 


30 


4.254 • 10 12 


3.895 • 10 12 


0.915 


40 


1.714 • 10 17 


1.630 • 10 17 


0.951 


50 


7.037 • 10 21 


6.825 • 10 21 


0.969 


60 


2.913 • 10 26 


2.857 • 10 26 


0.980 


70 


1.211 • 10 31 


1.196 • 10 31 


0.987 


80 


5.051 • 10 35 


5.007 • 10 35 


0.991 


90 


2.109 • 10 40 


2.096 • 10 40 


0.993 


100 


8.814 • 10 44 


8.774 • 10 44 


0.995 



4 Discussion 

4.1 A wider range of applicability 

The compensation approach for our counting problem leads to an exact expression for the generating 
function Q(x,y; z). A detailed exposition of the compensation approach can be found in [1, 3, 4], 
in which it has been shown to work for two-dimensional random walks on the lattice of the first 
quadrant that obey the following conditions: 

• Step size: Only transitions to neighboring states. 

• Forbidden steps: No transitions from interior states to the North, North-East, and East. 

• Homogeneity: The same transitions occur for all interior points, and similarly for all points 
on the horizontal boundary, and for all points on the vertical boundary. 

Although the theory has been developed for the bivariate transform of stationary distributions of 
recurrent two-dimensional random walks, this paper shows that is can also be applied to obtain a 
series expression for the trivariate function Q(x,y;z). The fact that Q(x,y;z) has one additional 
variable does not seem to matter much. We therefore expect that the compensation approach will 
work for walks in the quarter plane that obey the three conditions mentioned above. This will be 
a topic for future research. Another topic is to see whether the above conditions can be relaxed, 
and in fact, in the next section we present an example of a walk, not satisfying these properties, 
yet amenable to the compensation approach. 

4.2 Another example 

We consider the walk with S = {(— 1, 0), (— 1, — 1), (0, — 1), (1, — 1), (1, 0)} in the interior, Sh = 
{(-1,0), (1,0)} on the horizontal boundary, S v = {(0, 1), (0, -1), (1, -1), (1, 0)} on the vertical 
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boundary, and So = {(0, 1), (1,0)}; see the left display in Figure 5. This walk has a rather special 
behavior, as it can only move upwards on the vertical boundary. By simple enumeration, we get 
the following recursion relations: 



Qi,j/ Z = + + QiJ+1 + + Qi+l,j+li 

Qifi/z = qi-i,o + + %l + Qi+1,0 + 

Qo,j/z = ?o,j-i + Qoj+l + Qi,j + 
qofl/z = l/z + ^0,1 + + gi,o- 



i,j > 0, 
i > 0, 

i>o, 



(74) 
(75) 
(76) 
(77) 



Proposition 18. The unique solution to Equations (74) -(77) is given by q^j = coa Q (3 , with /3 
defined as the unique formal power series in z and solving 

p/z = l + /3 2 + /3 2 (l + /3) 2 , (78) 

raic (and even rational in x,y) and given 



a o — /3o(l + /3o) an d Co = Po/z. Therefore, the GF is 
by 



Q(x,y;z) 



oo 

E 

i,j,k=0 



co 



(1 - a x)(l - /3 y) 



(79) 



Proof. Let us first construct all pairs (a,/3) such that ca l j3 3 solves (74), (75) and (76). We obtain 
that a and /3 must satisfy 



a/3/z = p + /3 2 + a/3 2 + a 2 /3 + a 2 /3 2 , 
= a + a/3 2 + a 2 /3 + a 2 /3 2 . 



(80) 
(81) 



(Note that the reason why we obtain two and not three equations is that Equation (75) gives 
exactly the same information as (74).) Further, all pairs (a, /3) that satisfy both (80) and (81) are 
such that /3(1 + 0) = a. Substituting the latter into (80) gives (78). Denote by /3k, k G {0, . . . , 3}, 
the four roots of (78). Among the (3^, only one, say /3o, is a unique formal power series in z, see 
[20, Proposition 6.1.8]. Letting 

a = /3o(l + /3 ), (82) 

we conclude that for any cq, coa /3Q is a solution to (74), (75) and (76). With the particular choice 
co = 1/[1 — z(ao + ao/3o + /3o)], coa /3Q is also solution to (77). By uniqueness of the solution to 
the recursion relations (see Lemma 1), we conclude that qij = coa (3 J . Finally, by (78) and (82), 
we can simplify cq into (3q/z. □ 



Figure 5: Two walks with different boundary behavior 
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We now modify the boundary behavior of this walk and we show that this has severe conse- 
quences for the GF. Consider again the walk with S = {(— 1, 0), (— 1, — 1), (0, — 1), (1, — 1), (1, 0)} 
in the interior, but this time with Sy = {(0, — 1), (1, — 1), (1, 0)} on the vertical boundary, So = 
{(1,0)}, and the following rather special step set on the horizontal boundary: if the walk is in state 
(i, 0) with % > 0, the possible steps are {(—1, 0), (1, 0)} and a big step to (0, i); see the right display 
in Figure 5. We should note that a random walk, with similar unusual boundary behavior, has 
already been considered in [19]. The equations for q$j are then given by (74), (75), (77) and 



%,j/ z = 9j,o + 30,5+1 + Qi,j + 3 > 0. (83) 



Proposition 19. Let 



Po = 0, and define cq by 



2z(l + 1) 



Vco = E n+^'^i^ ' 1 - z{h + &+ 1 + Mk+l)] (85) 

^ (1 + Pt)--- (I + Pk) 

with 

/3 k+i = a k , a k+ i = g(Pk+i), c k+1 = c k ak + l , k ^ 0. (86) 

1 + Pk+i 

The unique solution to Equations (74), (75), (77) and (83) is given by 

oo 

$>4$- (87) 



Therefore, the GF is given by 



M,3 ~ 

k=0 



oo 

Q(x,y;z)=J2n % ^ ( 88 ) 

^ (1 - a k x){l - p k y) 



which has infinitely many poles in x and y, and is therefore nonholonomic as a trivariate function. 

Proof. Let us first notice that compared to Section 2, there is here a main novelty in the application 
of the compensation approach: it comes from the fact that Equation (74) for the interior and 
Equation (75) for the horizontal boundary are exactly similar. Accordingly, the compensation 
terms will appear by repeatedly considering Equation (76) for the vertical boundary only. 

We start from the product qij = ca 1 ^ . The pairs (ct,/3) that satisfy the inner equations (74) 
are characterized by a = g(f3) with g as in (84). (It is worth mentioning that a priori, we have two 
possibilities for g: one as in (84), and another with a sign + in front of the square root. We have 
to choose (84) because otherwise, the a k and (3 k in (86) would not be power series.) Substituting 
cofft into (83) yields 

ca j = cp j (l/z-a- fi-a/3). (89) 

This suggests to set /3q = 0, as otherwise (89) does not hold. (It is worth noting that this choice 
of (3q is less canonical than the choice we made in Lemma 2. However, if we are able to construct 
a solution to the recursion relations (74), (75), (77) and (83) with this particular choice, then by 
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uniqueness of the solution, it will be the only one, see Lemma 1.) Consider thus the product coagpQ, 
with /?o = and = g(/3o), for which (89) becomes 

c a = 0. (90) 

We then conclude that coOq/3q does not satisfy (83) and to correct for the error, we add a second 
product c\a\j3\, with (ai,/?i) on the curve a = g([3), in such a way that the term c\a\fi\ appears 
on the right-hand side of (90). This yields 

c a J Q = c 1 /3{(l/z-a- p-aP), (91) 

which gives 

1 OL\ 

Pi = a , ai=g(Pi), a = c — -3 5- = c , (92) 

l/z - ai - Px - oLiPx 1 + Pi 

where the last equality follows from (74). By adding the second product, we also obtain an extra 
term c\a\ on the left-hand side of (91). That is why we add a third product C2a\Pi' with (o^jp^) 
on the curve a = g(/3), and via similar reasoning, 

/3 2 = ai, a 2 =g(/3 2 ), c 2 = ci " 2 . (93) 

1 + P2 

Continuing this procedure means that we keep adding products with c^, a& and (3 k as in 

(86), so that (87) is a formal solution to Equations (74), (77) and (83). 
We now find cq. With (86), that we have proved above, we obtain 

Further, with (77) we have 



00 

co 



T - = - + co T ~(a k + p k + a k p k ). (95) 
* S Co Z to c ° 

Together with (86), we obtain (85). The a k and f3 k converge to the unique formal power series /3* 
solution to g(/3*) = Using (84), we obtain that 

/3, 3 + 2/32 + /3,(l-l/z) + l = 0. (96) 

The series /3* is such that /3*(0) = 0. It remains to show that q^j defined in (87) is a (well-defined) 
formal power series. There is here an interesting phenomenon: for the general compensation 
approach of Section 2, qij was a formal power series because a k = z 2k+1 6tk and fik = z 2k+2 (3k, 
where 2^ and (3k are formal power series (see Proposition 3). But here, we have otk = za k and 
fik = z(3k, where and /3& are formal power series not vanishing at (this follows from (84) and 
(86)). On the other hand, we have Ck = z k Ck, where Ck is a formal power series. Therefore, in the 
particular case of Proposition 19, qij is a formal power series thanks to the constants (see (86) 
and (87)). The nonholonomy of the trivariate function Q(x, y; z) follows by the same arguments as 
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that given below Proposition 8. 



□ 
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A Remaining proofs 

Let us start by expressing h, defined in (64), as an alternating sum. The reason why we wish to for- 
mulate h differently is twofold. First, this will enable us to show that p G [0.34499975, 0.34499976], 
see Corollary 22 — and in fact, in a similar way, we can approximate p up to any level of precision. 
Also, this is actually a key lemma for proving Lemma 13. 

Lemma 20. 



2011. 




(97) 



Proof. Equation (97) follows immediately from (49) and (64). 



□ 



As a preliminary result, we also need the following refinement of Proposition 10: 



Lemma 21. Assume that z G (0, l/y/8). Then both sequences {ak}k^o and {fik}k^o ar e positive 
and decreasing. Moreover, for all k ^ 




(98) 



Proof. It is enough to prove that for all 



t\ < 1 and z G (0,1/^8) 



(99) 



The proof of this inequality is similar to that of (53) and therefore omitted. □ 

Corollary 22. The series h(z) defined in (64) or (97) has a radius of convergence at least l/y/8. 
Its first positive zero p lies in [0.34499975,0.34499976]. 
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Proof. For z G (0,1/^8), the sequence {Tk}k^o denned in Lemma 20 is positive and decreasing, 
see Lemma 21. In particular, denoting 

A p = ]T(-l) fc T fc (100) 

k=0 

and using (97), for all p ^ and all z G (0, l/v^) we have 

1 + 2z(-l - a + A 2p+1 ) < h(z) < 1 + 2z(-l - a + A 2p ). (101) 

Applying the last inequalities to p = 4 and noting that the right-hand side (resp. left-hand side) 
evaluated at 0.34499976 (resp. 0.34499975) is negative (resp. positive) concludes the proof. □ 

Proof of Lemma 13. First, using the explicit expression of and (3^, and employing calculus soft- 
ware, we obtain that the algebraic function 1 + 2z(— 1 — ao + A 5 ) has only one zero within the circle 
of radius l/\/8, and that 

inf |l + 2z(-l -a + A 5 )| > 10~ 2 . (102) 

|*|=i/V5 

Note that the lower bound 10~ 2 in (102) is almost optimal; it is rather small because like p, the 
unique zero of 1 + 2z(— 1 — ao + A 5 ) within the disc of radius l/\/8 is very close to 1/V8. 

By Rouche's theorem (see Theorem 11), applied to the circle with radius l/\/8, it is now enough 
to prove that 

sup |[l + 2z(-l -a + A 5 )] - [l + 2z(-l -q + A°°)]| < 10~ 2 . (103) 

\z\=i/VE 



For this we write 

-2p+5 



|A°° — A p | 



E < E mi < T iZ i§^ < vv?> +3 , <™> 



where the second upper bound is obtained from the inequality jT^j ^ l/\/2 2fc+3 , see (97) and 
Lemma 21. By (104) we then obtain 

|[1 + 2z (-l - a + A p )] - [1 + 2z(-l - a + A°°)]| < 2|*|— (105) 

If |z| < l/v8, the last quantity is bounded from above by 1/2 P+2 . For p = 5, we obviously get 
1/2 P+2 < 10~ 2 , and (103) is proven. □ 

Lemma 23. 

C 0i0 = 0.0531 • [1 +O(10~ 3 )]. (106) 
Proof. As in the proof of Proposition 15, we can write 

C0,0 = —rZ < -^T, f yj" ) (107) 

Pi 1 - P x o,o(P)\ 
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so that the main difficulty lies in approximating x' 00 (p). For this, we use (49) and (97) to write 

oo 

£&,„(*) = -2a' (z) + J>l) fc T^). (108) 

fc=0 

While it is easy to control the series X]fcLo( — ^-) k ^k( z ), because a k (z) and /3 k (z), and hence T k (z), 
decrease exponentially fast to as k — > oo, see Proposition 3 and Lemma 21, it is not obvious 
how we should deal with Y^/k=o(~ -O^feC 2 )' where terms like a' k (z) and (3' k (z) appear. We next 
show that a' k (z) and (3' k (z) actually also decrease exponentially fast to as k — > oo, at least for 
z G [1/4,0.35]. Note that this assumption on z is not restrictive, since p belongs to the interval 
[1/4,0.35] by Corollary 22. 

Consider the sequence {7fc(^)}fc>o defined by 70(2) = ao(z) and, for k ^ 0, by 7fc + i(z) = 
f(lk{z)), with cto(z) and / as defined in (20) and (29), respectively. Note that (29) gives 72fc(-z) = 
a k (z) and 72fc+i(-z) = h{z). The sequence {^' k {z)}k^o satisfies the recurrence relation 

Yfc+iC*) = i' h {z)d t f(n k {z)) + d z f{ lk {z)). (109) 

Below we study the coefficients <9t/(7fc(z)) and d z f{^ k {z)) of (109). By using expression (29) of /, 
we easily obtain that for all z G [0, l/v^] and all t G [0, ao(.z)] 

\d t f{lk{z))\ < 4V2/9. (110) 

In addition, d z f(t) = f(t)/[zyjl - 4z 2 (l + 1 2 )], in such a way that \d z f{t)\ ^ \2/z\ ■ \f(t)\ for all 
z G [0, l/\/8] and all t G [0, ao(z)]. Using then Lemma 21, we obtain 

\d z f( 7k (z))\ < \2/z\ ■ \f(lk(z))\ = \2/z\ ■ \jk+i(z)\ ^ \2/z\ ■ 1/V2 k+1 < (111) 
where the last inequality follows from the assumption z G [1/4,0.35]. From (109)-(111) we get 

K+i(*)K'HYfc+i(*)l+& J ^ = 4^2/9, £ k = 8/V2 k+1 . (112) 

We deduce that 

k 

|Y fc+1 (*)| < r? fc+1 • | 7o (*)| + < • |V (z)| + 72/v^ fc+1 < 100/v^ + \ (113) 

p=0 

where the last inequality follows from sup 2g r .25,0.35] ^oC- 2 )! ^ ^ an< ^ V ^ Vv^- With Lemma 21 
and (113), we obtain 

|2£(z)| < 200/v / 2"' +2 , (114) 

and thanks to (108) we can obtain approximations of x' 00 (p) up to any level of precision. Lemma 23 
follows. □ 
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